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PREFACE 


This  report  describes  a computer  model  of  an  electric  trolley 
bus  developed  at  the  Transportation  Systems  Center  to  simulate  the 
power/propulsion  characteristics  of  an  urban  trolley  bus.  This 
model  comprises  the  fourth  in  a series  of  models  being  developed 
at  TSC  to  simulate  different  types  of  bus  propulsion  systems. 
Earlier  reports  describing  these  models  include:  1)  "Flywheel/ 
Diesel  Hybrid  Power  Drive:  Urban  Bus  Simulation,"  DOT-TSC-78-10, 
UMTA-MA- 06- 0044-78-1 , Final  Report,  May  1978;  2)  "Diesel  Bus 
Performance  Simulation  Program,"  DOT-TSC-UMTA-79-16,  UMTA-MA- 06- 
0044-79-1,  Final  Report,  April  1979;  and  3)  "Simulation  of  an 
Urban  Battery  Bus,"  DOT-TSC-UMTA-79-15,  UMTA-MA- 06- 0093-7 9- 1 , 

Final  Report,  July  1979,  These  simulation  models  provide  the 
capability  for  rapid  evaluation  of  bus  performance  through  the 
comparison  of  simulated  bus  performance  with  data  obtained  from 
engineering  tests.  The  work  conducted  in  this  area  is  sponsored 
by  the  Urban  Mass  Transportation  Administration  and  is  part  of  a 
larger  effort  concerned  with  the  demonstration  and  test  evaluation 
of  advanced  bus  propulsion  concepts. 
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EXECUTIVE  SUMMARY 


This  report  describes  a computer  model  developed  at  the 
Transportation  Systems  Center  to  simulate  power/propulsion  char- 
acteristics of  an  urban  trolley  bus.  This  model  comprises  the 
fourth  in  a series  of  models  being  developed  at  TSC  to  simulate 
different  types  of  bus  propulsion  systems.  The  work  conducted  in 
this  area  is  sponsored  by  the  Urban  Mass  Transportation  Administra- 
tion as  part  of  a larger  effort  concerned  with  the  demonstration 
and  test  evaluation  of  advanced  bus  propulsion  concepts.  The 
computer  models  developed  by  TSC  provide  a means  of  rapidly 
evaluating  bus  performance  and  of  comparing  simulated  bus  perfor- 
mance with  results  obtained  from  engineering  tests.  They  also 
provide  a means  of  developing  comparative  information  on  the 
performance  characteristics  and  component  sizing  requirements  of 
various  bus  propulsion  systems.  The  results  of  the  computer  an- 
alyses will  be  incorporated  into  the  Annual  Technological  Assess- 
ment Report  for  future  reference  and  evaluation. 

This  report  includes  a description  of  the  application  of  the 
trolley  bus  model  to  a typical  urban  bus  drive  cycle.  The  re- 
quested mission  profile  is  given  as  input  to  the  bus  model  while 
the  actual  mission  profile  determined  by  the  limitation  of  the  bus 
propulsion  system  appears  as  output  data.  The  output  data  file 
also  lists  the  electrical  drive  characteristics,  including  com- 
ponent power  losses,  at  successive  intervals  of  drive  time. 
Applications  of  the  trolley  bus  model  anticipated  for  the  future 
include  the  evaluation  of  energy  consumption,  the  impact  of  com- 
ponent sizing  on  bus  operation,  and  the  effect  of  mission  profile 
on  the  ultimate  efficiency  of  the  trolley  bus.  Studies  of  this 
type  provide  a useful  data  base  for  assessing  the  capabilities  of 
the  bus  propulsion  system  in  urban  transit  applications. 
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1,  INTRODUCTION 


This  report  describes  a computer  model  of  the  propulsion 
system  of  an  electric  trolley  bus  and  the  application  of  this 
model  for  simulating  trolley  bus  performance.  The  computer  model 
has  been  implemented  on  the  DECsystem-10  computer  located  at  the 
Transportation  Systems  Center  (TSC)  in  Cambridge,  Massachusetts. 
Written  in  FORTRAN-10  source  language,  the  program  could  be  run  on 
most  large-scale  computer  systems. 

The  program  models  an  electric-powered  passenger  vehicle 
traversing  a defined  mission  profile  of  acceleration,  cruising 
speed,  roadway  grade,  and  encountered  headwind.  Sufficient 
flexibility  has  been  introduced  into  the  system  (Model)  to  allow 
for  arbitrary  mission  profiles  and  variations  in  the  system  com- 
ponent specifications.  The  program  has  been  designed  to  evaluate 
the  power  consumption  of  a trolley  bus  under  varied  conditions  of 
operation,  including  the  use  of  regeneration  for  optimum  power 
utilization . 

The  computer  simulation  of  the  electric  trolley  bus  has 
several  objectives.  First,  it  offers  a convenient  means  of  assess- 
ing the  performance  capabilities  of  an  electric  trolley  bus  and  of 
comparing  its  performance  writh  that  of  other  electric  and  non- 
electric vehicles.  Such  comparisons  are  useful  for  evaluating  the 
relative  merits  of  different  vehicle  propulsion  systems.  The 
conditions  under  which  one  system  can  be  optimized  relative  to 
another  system  are  easily  explored  via  computer  models.  Second, 
the  computer  simulation  affords  a means  of  examining  the  impact 
of  a particular  speed-time  profile  on  trolley  bus  performance. 

The  choice  of  one  transit  profile  over  another  can  be  made  using 
the  computer  program  to  estimate  energy  losses  under  varying  pro- 
file conditions.  The  effect  of  power  regeneration  on  overall 
vehicle  efficiency  as  well  as  the  best  choice  of  power  control 
equipment  for  optimum  vehicle  performance  can  also  be  easily 
assessed  using  computer  simulation  techniques. 
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The  purpose  of  this  report  is  to  document  the  electric 
trolley  bus  program  and  to  provide  the  reader  with  the  necessary 
information  to  operate  the  program  and  to  interpret  the  results 
obtained  from  its  use.  Detailed  procedures  are  given  in  the  form 
of  a user's  guide  for  executing  the  program  on  a DECsystem-10 
computer.  These  are  summarized  in  Section  3.2  of  the  report. 
Section  5 contains  FORTRAN  source  listings  of  the  main  program 
and  all  subroutines  as  well  as  sample  output  data.  As  a guide 
for  understanding  the  model  and  interpreting  the  computer  results, 
the  modeling  theory  and  equations  are  detailed  in  Section  4. 

1.1  MODULAR  PROGRAM  FORMAT 

The  computer  model  utilizes  a modular  program  format  in 
which  functions  common  to  a given  system  component  are  contained 
in  a separate  subroutine.  For  example,  all  operations  involving 
the  generation  of  tractive  effort  and  flow  of  power  to  and  from 
the  traction  machine  are  delegated  to  the  subroutine  MOTOR.  The 
advantage  of  this  format  is  that  it  facilitates  the  interchange  of 
alternative  types  of  system  components  without  impairing  the  basic 
integrity  of  the  overall  program.  This  is  particularly  useful 
when  one  wishes  to  explore  the  effect  of  changing  the  design 
parameters  of  a system  component  or  to  examine  the  effect  of  a 
particular  mode  of  operation  on  the  propulsion  performance.  This 
format  has  been  utilized  in  previous  modeling  efforts  (e.g., 
battery  bus,  AGRT  systems)  with  favorable  results. 

Figure  1 shows  a block  diagram  of  the  modular  format  as 
applied  to  the  electric  trolley  bus.  The  separate  subroutines 
which  perform  specific  functions  are  PROFILE,  PCU,  and  MOTOR. 

The  MAIN  program  provides  the  calling  statements  and  sequence  logic 
necessary  to  maintain  program  continuity. 

1.2  TROLLEY  BUS:  RESEARCH  CONTEXT 

The  renewed  interest  in  electric  trolley  vehicles  has  been 
stimulated  by  several  factors.  These  rubber-tired  transit  vehicles 
are  quiet,  waste  little  energy  at  idle,  and  emit  no  exhaust  fumes. 
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FIGURE  1.  SKETCH  OF  MODULAR  PROGRAM  FORMAT 
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They  are,  therefore,  well  adapted  to  the  great  majority  of  urban 
driving  environments.  Because  of  the  possibility  of  built-in 
storage  devices,  they  have  the  flexibility  of  operating  for  cer- 
tain distances  off-wire. 

Electric  vehicles  possess  better  acceleration  characteris- 
tics, even  under  load,  than  do  internal  combustion  vehicles. 
Consequently,  they  are  capable  of  doing  more  work  in  a given  time 
period  compared  writh  other  vehicles.  Furthermore,  because  of  the 
scarcity  and  increasing  costs  of  petroleum-based  fuels,  alternative 
forms  of  energy  propulsion  systems  are  becoming  increasingly  more 
attractive.  Since  the  supply  of  fossil  fuels  far  exceeds  that  of 
petroleum-based  fuels,  vehicles  powered  by  electricity  generated 
from  fossil  fuels  are  becoming  more  promising  for  the  future. 

Programs  concerned  with  the  development  of  electric-powered 
vehicles  are  now  under  way  in  a number  of  transportation- or iented 
companies  in  the  United  States.  These  include  among  others, 

General  Electric,  Boeing,  Garrett  Airesearch,  Vought , and 
Westinghouse . These  programs  involve  the  development  of  electric- 
drive  vehicles  ranging  in  size  from  small  commuter  cars  to  large 
mass  transit  buses.  In  addition,  studies  are  now  in  progress  to 
determine  the  practicality  of  battery-powered  buses  and  supplemen- 
tary energy  storage  schemes,  including  flywheels. 
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2,  BASELINE  SYSTEM:  COMPONENTS  AND  MODES  OF  OPERATION 


2.1  BASELINE  SYSTEM 

The  baseline  system  comprises  a 24,400  lb  trolley  coach 
(fully  loaded)  driven  by  a GE  1213,  123  hp  DC  traction  motor.  A 
600  volt  overhead  trolley  supplies  power.  An  MRC  (multiple  relay) 
controller  provides  control  of  traction  motor  speed  by  means  of 
accelerator  and  brake  pedals  on  the  floor.  The  propulsion  system 
is  capable  of  accelerating  the  coach  at  3 mph/sec  up  to  a speed  of 
17  mph  and  at  reduced  accelerations  up  to  the  peak  vehicle  speed 
of  49  mph.  Vehicle  deceleration  is  achieved  through  either 
regenerative  or  dynamic  frictional  braking. 


2.2  CONFIGURATION  OF  MAJOR  COMPONENTS 

The  major  components  of  the  Trolley-Bus  Simulation  include 
the  vehicle  coach,  DC  compound  motor,  MRC  (cam)  controller,  and 
power  source  (catenary) . A block  diagram  showing  the  system 
components  is  given  in  Figure  2.  The  important  quantities 
associated  with  each  component  are  indicated  below  the  diagram. 
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Control 

Resistance 


* armature 
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armature 
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FIGURE  2.  BLOCK  DIAGRAM  OF  TROLLEY  BUS  SYSTEM 


2.2.1  Trolley  Coach 

The  trolley  bus  modeled  in  the  computer  program  is  a rubber- 
tired  mass  transit  coach  with  a seating  capacity  of  15  passengers. 
Peak  speed  under  normal  operating  conditions  is  45  mph.  Table  1 
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summarizes  the  coach  parameters  used  in  the  computer  model. 


TABLE  1.  TROLLEY  COACH  SPECIFICATIONS 


Trolley  Coach  Specifications 

Curb  Weight  (WrTCRB) 

= 22,000 

lb 

Gross  Weight  (15  passengers) 

= 24,400 

lb 

Frontal  Area  (AF) 

= 69.5  ft 

2 

Aerodynamic  Drag  Coeff.  (CD) 

= 0.66 

Rolling  Drag  Coeff.  (CR) 

= .0075 

Tire  Factor  (TF) 

= 8.003 

Gear  Ratio  (RATIO) 

= 11.59  reduction 

Wheel  Diameter 

= 4 2 in 

2.2.2  Vehicle  Mission  Profile 

The  "mission"  profile  describes  vehicle  speed  as  a function  of 
time.  It  can  either  be  defined  by  a specified  vehicle  speed-time 
dependence  or  derived  from  input  data  on  vehicle  speed  and  accelera- 
tion (and  deceleration).  The  latter  approach  is  the  one  used  in 
this  simulation  study. 

Figure  3 shows  one  complete  driving  cycle  of  an  arbitrary  mis- 
sion together  with  certain  key  mission  profile  parameters.  The 
vehicle  traverses  a simulated  bus  route,  which  is  defined  by  the 
number  of  stops  per  mile  (NS)  and  the  route  length  (SR) . Multiple 
transits  of  the  same  route  can  also  be  run,  with  an  optional  wait- 
ing period  added  between  routes  (TAUMU) . During  this  waiting  period 
and  whenever  the  vehicle  is  at  rest,  only  the  auxiliary  loads 
(including  environmental  control  systems)  use  power.  The  mission 
profile  is  further  defined  by  the  cruise  speed  between  stops  (VC) , 
the  required  acceleration  (ACCM)  and  deceleration  (DECL) , the 
encountered  head  wind  (HW) , and  the  roadway  grade  (GRADE) . The 
stops  are  assumed  to  be  equally  spaced,  and  equal  dwell  time  (DWELL) 
is  spent  at  each  stop  for  passenger  boarding.  This  route  model 
generates  a repeating  series  of  driving  cycles  between  each  stop 
as  the  vehicle  traverses  the  mission. 
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FIGURE  3.  SPEED-TIME  PROFILE  FOR  ONE  CYCLE  OF  AN 
ARBITRARY  MISSION 


The  mission  profile  is  "specified"  at  the  start  of  each  com- 
puter run  by  selecting  one  of  three  programmed  driving  cycles, 
called  cycles  A,  B,  and  C.  The  values  of  the  different  program 
constants  associated  with  the  A,  B,  and  C drive  cycles  are  given  in 
Table  2.  These  drive  cycles  describe  the  vehicle  speed  and 
acceleration  as  "requested"  by  the  bus  operator.  Actual  vehicle 
speed  and  acceleration  are  determined  by  the  available  motor  out- 
put power  as  well  as  by  limitations  in  the  control  characteristics 
of  the  power  control  unit  (PCU). 

The  vehicle  starts  from  rest  in  a jerk-limited  mode,  (See 
Figure  3.)  As  the  acceleration  (ACC)  increases  to  the  required 
rate  (ACCM) , the  jerk  (AJERK)  (which  is  the  rate  of  change  of 
acceleration)  is  constrained  by  a specified  limit  (AJERKR) . As 
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the  vehicle  approaches  cruise  speed  (VC),  a j erk- limited  decrease 
of  acceleration  is  programmed.  At  the  proper  speed,  the  vehicle 
acceleration  is  reduced  to  zero  at  the  jerk-limited  rate.  During 
cruise,  the  controller  (CAM)  seeks  to  maintain  constant  cruise 
speed  by  alternately  applying  positive  and  negative  tractive 
efforts.  During  deceleration,  the  vehicle  undergoes  a programmed 
deceleration  profile  in  a reverse  sequence  to  the  acceleration 
profile.  The  vehicle  is  brought  to  rest  with  a specified 
deceleration  (DECL) . 


TABLE  2.  DRIVE  CYCLE  PARAMETER  VALUES 


Driving  Cycle 

Data 

DESCRIPTION 

FORTRAN 

SYMBOL 

CYCLE 

A 

CYCLE 

B 

CYCLE 

C 

UNITS 

Thrust  limited  acceleration 

ACCM 

3.5 

3.5 

3.0 

mph/sec 

Deceleration  rate 

DECL 

1.5 

3.5 

3.5 

mph/sec 

Cruise  speed 

VC 

20.0 

25.0 

25.0 

mph 

2 

Maximum  allowed  jerk 

AJERKR 

3.5 

3 . 5 

3.5 

mph/sec 

Number  of  stops  per  mile 

NS 

5 

8 

10 

Distance  between  stops 

SS 

0.20 

0 . 125 

0.10 

mi 

Dwell  time  at  each  stop 

DWELL 

12.0 

16.0 

8.0 

sec 

Additional  dwell  at  end 
of  route 

TAUMU 

0.0 

0.0 

30.0 

sec 

Number  of  passengers 

NP 

20 

10 

15 

Route  length 

SR 

6.0 

6.0 

6.0 

mi 

Number  of  transits  of  route 

NT 

4 

4 

4 

Armature  current  setpoint 

SET 

200.0 

250.0 

225.0 

amp 

2.2.3  DC  Traction  Motor 


A GE  1213^  compound  wound,  150  hp , 600  volt  DC  traction  motor 
is  used  in  the  trolley  bus  simulation.  The  electrical  specifica- 
tions of  the  motor  are  given  in  Table  3. 
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TABLE  3.  DC  MOTOR  (GE  1213)®  ELECTRICAL  SPECIFICATIONS 


Motor  Electrical  Specifications 

Type:  GE  1213  Compound  Wound,  Commutating  Pole 

No.  of  Poles  = 4 
Ratings : 

Armature  voltage  = 600  vac 
Armature  voltage  (AVE)  = 550  vac 
Armature  current  = 205  A 
Shunt  field  current  = 2.5  A 
Output  Power  (Shaft)  = 150  hp 
Max.  speed  = 4500  rpm 

Winding  Data: 

No.  of  shunt  field  turns  = 1175 
No.  of  series  field  turns  = 12 
Armature  resistance  = .05140  min., 

Shunt  field  resistance  = 80.50  min., 

Series  field  resistance  = .01250  min, 
Commutating  resistance  = .0220  min., 


05540  max. 
85.50  max . 

, .01350  max, 

.0300  max . 


Additional  information  pertaining  to  the  physical  characteris- 
tics of  the  GE  1213  motor  is  given  in  Table  4. 

TABLE  4.  MOTOR  (GE  1213)  MECHANICAL  SPECIFICATIONS 


Weight : 

Motor  (total) 

1429 

lb 

Armature 

395 

lb 

Armature  Moment  of  Inertia 

35 

lb-ft2 

Max.  Dynamic  Braking  Capability 

2400 

lb 
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Figure  4 shows  the  measured  tractive  effort  (at  wheels)  and 
vehicle  speed  as  a function  of  the  DC  armature  current.  Tractive 
force  increases  almost  linearly  with  armature  current  for  a given 
shunt  field  excitation.  At  vehicle  speeds  above  14  mph,  it  is 
necessary  to  reduce  the  shunt  field  excitation  (f ield-weaken) , as 
indicated  in  the  figure. 

Regenerative  braking  can  be  accomplished  at  higher  speeds  by 
increasing  the  motor  field  excitation  sufficiently  to  cause  the 
back  EMF  to  exceed  motor  terminal  voltage.  The  resulting  braking 
effort  (at  the  wheels)  realized  for  different  shunt  field  excita- 
tions is  shown  in  Figure  5.  Dynamic  (rheostatic)  braking  can  be 
accomplished  over  the  full  range  of  vehicle  speeds  by  disconnecting 
the  input  power  to  the  motor  and  connecting  a suitable  resistor 
grid  across  the  armature  terminals.  Figure  5 show's  the  braking 
effort  developed  with  a 1.5  ohm  braking  resistor  and  three  different 
field  excitations. 

Magnetic  saturation  plays  an  important  part  in  determining 
motor  performance.  Details  regarding  the  manner  in  which  magnetic 
saturation  is  included  in  the  motor  simulation  are  discussed  in 
Section  4.2. 

2.2.4  Power  Controller  (PCU) 

Power  to  the  DC  motor  is  controlled  by  a type  MRCl  controller 
which  the  vehicle  operator  activates  by  pushing  down  on  the 
accelerator  pedal.  As  the  accelerator  pedal  is  depressed,  the 
controller  cam  assumes  1 of  14  different  step  positions.  At  each 
step  position,  a predetermined  resistance  is  inserted  into  the 
armature  and  shunt  field  circuits.  Table  5 lists  the  different 
'cam'  resistances  and  percent  normal  shunt  field  current  for  the 
various  step  positions.  Inserting  resistance  into  the  armature 
circuit  controls  the  motor  terminal  voltage,  wrhile  inserting 
resistance  into  the  shunt  field  circuit  gives  field  (flux)  control. 
Cam  positions  1 through  10  give  combined  terminal  voltage  and  field 
control;  cam  positions  10  through  14  yield  field  control. 
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CHARACTERISTIC  CURVES  ON  S50  VOLTS 


CALCULATED  FROM  TESTS 
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FIGURE  5.  GE  1215^  TRACTION  MOTOR  BRAKING  CHARACTERISTICS 


TABLE  5.  MRC  CONTROLLER  RESISTANCE,  CURRENT  CHARACTERISTICS 


Cam 

Position 

Accelerating 

Percent  of 
Total  Series 
Resistance  in 
Series  with 
Motor 

Percent  Nor. 
Shunt  Field 
Current 

Approx.  Bal. 
Speed  MPH 
550  Volts 

Open  Circuit 

0 

1 

100 

0 

2 

100 

26 

3 

68 

33 

4 

50 

100 

5 

39 

100 

6 

27 

100 

7 

20 

100 

8 

12 

100 

9 

5 

100 

10 

0 

100 

17 

11 

0 

72 

20 

12 

0 

46 

24 

13 

0 

24 

30 

14 

0 

0 

40 

Regen.  Brake 

13  to  10  inc . 

0 

24  to  1001 

40  to  17 

Dynamic  Braking 

1 

52 

73 

2 

52 

100 

3 

52 

112 

4* 

52 

112 

^Gradual  increase  from  step  3 to  4 (full) . 
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The  MRC  controller  functions  equally  in  both  motoring  and 
regenerative/dynamic  braking  modes.  In  the  regenerative  braking 
mode,  the  motor  back  EMF  must  exceed  the  terminal  voltage.  For 
a terminal  voltage  of  550  volts,  this  requires  that  the  vehicle 
speed  be  greater  than  18  mph.  (See  Figure  5.)  To  produce  regen- 
erative braking,  the  vehicle  operator  releases  pressure  on  the 
power  (accelerator)  pedal,  allowing  the  cam  to  assume  a lower 
position  within  the  10-14  step  range  of  the  cam.  Figure  5 shows 
the  braking  tractive  effort  realized  for  three  values  of  shunt  field 
current,  while  Figure  4 gives  the  corresponding  armature  current 
for  different  tractive  efforts.  Since  the  tractive  force  is 
negative  in  this  case,  armature  current  is  likewise  negative. 

For  dynamic  (rheostatic)  braking,  the  vehicle  operator  shifts 

control  from  the  accelerator  pedal  to  the  brake  pedal.  Pushing 

dowm  on  the  brake  pedal  causes  the  motor  to  be  disconnected  from 

the  line  and  reconnected  to  a brake  resistor.  The  DC  machine  then 

functions  as  a generator,  converting  the  vehicle's  kinetic  energy 
2 

into  I R heat  loss  in  the  brake  resistor.  To  supplement  the 
braking  effort,  particularly  at  low  vehicle  speeds,  friction 
braking  is  added.  The  way  in  which  the  power  and  brake  pedals  are 
related  to  vehicle  tractive  effort  is  illustrated  in  Figure  6. 

2.3  MODES  OF  OPERATION 

Block  diagrams  of  power  flow  in  motoring  and  braking  modes  of 
operation  are  shown  in  Figure  7.  In  the  motoring  modes,  the  power 
flows  from  the  trolley  wire  through  the  power  conditioning  unit 
(PCU)  to  the  traction  motor.  The  output  motor  power  then  flows 
through  the  gearbox  and  drivetrain  to  the  wheels.  In  the 
regenerative  braking  mode,  the  power  flow  is  reversed  in  the  main 
path.  The  controller  maintains  cruise  speed  by  alternating 
between  motoring  and  regenerative  braking,  with  the  net  accelerat- 
ing thrust  averaging  out  to  zero.  Service  braking  is  accomplished 
by  a blend  of  dynamic  and  friction  brakes.  Three  rates  of 
deceleration  are  available  from  the  dynamic  braking  system,  corres- 
ponding to  three  values  of  shunt  field  resistance.  Any  decelera- 
tion rate  is  possible,  however,  by  assuming  an  ideal  blending  of 


14 


FIGURE  6. 


SKETCH  INDICATING  POWER  AND  BRAKE  PEDAL  STEP  POSITIONS 
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FIGURE  7.  BLOCK  DIAGRAM  SHOWING  POWER  FLOW  IN 
MOTORING  AND  BRAKING  MODES 
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friction  brakes.  During  dynamic  braking,  the  PCU  controls  only 
the  shunt  field  and  auxiliary  power,  which  is  supplied  from  the 
line . 

During  acceleration,  the  armature  current  falls  as  the  motor 
back  EMF  builds  up.  To  maintain  a constant  acceleration,  the 
vehicle  operator  "selects"  a higher  cam  position,  which  causes  an 
abrupt  increase  in  armature  current.  From  standstill  to  the  final 
cruise  velocity  condition,  the  motor  armature  current  undergoes 
successive  abrupt  transitions  corresponding  to  the  increasing  cam 
positions.  The  manner  in  which  the  armature  current  varies  in 
response  to  "requested"  acceleration  is  shown  in  Figure  8 for 
the  drive  cycle  C profile.  The  cam  position  step  number  is  in- 
dicated in  the  figure  for  each  interval  in  the  drive  cycle. 

During  cruise,  the  armature  current  is  reduced  in  response  to 
the  lower  tractive  effort  required  to  overcome  the  vehicle  drag 
forces.  Due  to  the  coarseness  of  the  controller  action,  the  trolley 
bus  model  predicts  alternately  positive  and  negative  tractive 
efforts  as  necessary  to  maintain  a constant  cruise  velocity.  This 
results  in  alternately  positive  and  negative  armature  currents, 
as  shown  in  Figure  8. 

The  deceleration  mode  uses  regenerative,  dynamic,  and  fric- 
tional braking  to  reduce  the  vehicle  speed.  In  the  regenerative 
and  dynamic  braking  modes,  motor  armature  current  flows  in  the 
negative  (reverse)  direction  as  required  to  develop  braking  torque. 
For  the  case  illustrated  in  Figure  8,  in  which  dynamic  and  fric- 
tional braking  are  jointly  utilized  for  decelerating  the  vehicle, 
the  armature  current  shown  is  the  current  flowing  through  the 
braking  resistor. 

The  choice  of  the  particular  braking  mode  used,  whether 
regenerative,  dynamic  and/or  frictional,  is  dictated  by  the  amount 
of  braking  effort  required  for  deceleration.  Initially,  regenera- 
tive braking  is  tried  by  the  vehicle  operator.  If  additional 
braking  effort  is  required,  the  operator  removes  his  foot  from 
the  motoring  pedal  and  pushes  down  on  the  brake  pedal.  This  results 
in  a combination  of  dynamic  and  frictional  braking  sufficient  to 
meet  the  required  braking  demands. 
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MOTOR  ARMATURE  CURRENT 


FIGURE  8.  MOTOR  ARMATURE  CURRENT  FOR  DIFFERENT 
OPERATING  MODES 
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3,  TROLLEY  BUS  PERFORMANCE  PROGRAM 


The  Trolley  Bus  Performance  Program  consists  of  a main  pro- 
gram and  three  subroutines.  The  main  program  controls  the  execution 
and  calls  the  subroutines  in  sequence.  Each  subroutine  contains 
modeling  equations  which  determine  the  operating  characteristics 
of  a part  of  the  propulsion  system. 

The  subroutine  PROF  determines  the  mission  profile  conditions 
at  each  program  iteration.  These  conditions  include  acceleration, 
speed  and  position  of  the  vehicle,  roadway  grade,  and  encountered 
headwind.  PROF  also  computes  the  distance  from  the  vehicle  to  the 
next  stop  and  determines  when  to  begin  braking.  Braking  is  done 
at  a constant,  but  j erk -1 imited , deceleration  rate  that  is  specified 
by  the  driving  cycle. 

The  subroutine  MOTOR  solves  the  motor  modeling  equations  and 
determines  the  motor  speed,  armature  and  field  current,  generated 
emf  and  terminal  voltage,  power  losses,  and  power  input  to  the 
propulsion  system  from  the  trolley  wires. 

The  subroutine  PCU  computes  the  power  lost  in  the  dc/dc 
power  conditioning  unit  and  the  motor  field  controller,  and  also 
determines  the  power  supplied  to  auxiliary  systems. 

The  computer  program  calculates  and  prints  out  numerous  system 
parameters  at  each  time  increment.  Profile  conditions,  tractive 
force,  motor  parameters,  and  vehicle  power  and  energy  are  printed 
in  the  data  file  FOR05.DAT.  At  the  end  of  the  run,  a run  summary 
is  also  printed  in  FOR03.DAT  which  contains  the  distance,  time, 
and  average  speed  of  the  mission,  as  well  as  the  total  energy 
consumed  and  the  energy  consumed  per  mile. 

A plotting  program  has  also  been  developed  to  graphically 
display  variable  system  parameters  in  six-channel  stripchart 
fashion  on  an  off-line  CALCOMP  pen  plotter.  The  plotting  capabil- 
ity is  further  discussed  in  Section  3.3  of  this  report. 
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3.1  PROGRAM  FLOW  AND  SUBROUTINE  CALLING  SEQUENCE 


A simplified  flowchart  of  the  computer  program  is  shown  in 
Figure  9.  Complete  FORTRAN  source  listings  are  given  in  Section 
4.0. 

The  subroutine  PROF  is  called  first  to  initialize  parameters 
within  that  subroutine.  On  all  subsequent  calls  to  PROF,  the 
initialization  section  is  skipped.  Page  and  column  headings  are 
then  written  in  the  output  data  file  FOR03.DAT,  and  the  program 
enters  its  main  iteration  loop. 

At  each  time  increment,  the  subroutine  PROF  is  called  to 
determine  the  mission  profile  conditions  of  acceleration,  speed 
and  position  of  the  vehicle,  roadway  grade,  and  encountered  head- 
wind. PROF  also  computes  the  total  drag  force  on  the  vehicle  and 
determines  when  to  begin  braking  so  that  the  vehicle  comes  to 
rest  at  the  correct  position.  The  subroutine  MOTOR  is  then  called 
to  compute  the  motor's  performance  characteristics  and  power  input. 
The  normal  acceleration  characteristic  is  obtained  by  maintaining 
the  motor  armature  current  at  or  above  a specified  setpoint.  This 
is  done  by  removing  series  resistors  from  the  armature  circuit  in 
steps.  There  are  14  such  steps  corresponding  to  positions  on  a 
rheostatic  cam  controller.  The  subroutine  PCU  is  then  called  to 
compute  the  PCU  power  losses  and  auxiliary  load.  The  total  power 
input  to  the  vehicle  and  the  total  energy  consumed  since  the  start 
of  the  run  are  then  computed. 

The  program  then  writes  one  'line  of  data  in  the  output  file 
FOR03.DAT.  If  plots  are  desired,  one  line  of  data  is  also  written 
in  the  file  FOR22.DAT.  (This  file  is  required  as  input  to  the 
plotting  program  and  contains  the  values  of  six  time-varying 
system  parameters  that  are  to  be  displayed.)  The  time  clock  is 
then  incremented  if  the  end  of  the  cycle  has  not  been  reached,  and 
the  program  loops  back  to  PROF  for  the  next  iteration. 

A specified  number  of  driving  cycles  are  made  per  mile;  and 
at  the  end  of  each  cycle  the  stop  index  is  incremented,  page  and 
column  headings  written,  and  the  next  cycle  is  begun.  At  the  end 
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FIGURE  9.  PROGRAM  FLOWCHART 
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of  a specified  number  of  routes  the  mission  is  complete  and  the 
run  summary  calculations  are  performed  and  the  results  written  in 
FOR03.DAT. 

The  plotting  program  CHART. FOR  can  then  be  executed  to  dis- 
play any  six  time-varying  system  parameters  on  an  off-line  CALCOMP 
pen  plotter. 

3.2  EXECUTING  THE  PROGRAM 

The  main  program  and  all  subroutines  are  contained  in  separate 
source  files.  By  not  combining  all  subroutines  into  one  file, 
compilation  time  is  minimized  in  case  changes  must  be  made  to  only 
a few  files.  The  main  program  and  the  three  subroutines  must  be 
compiled  and  loaded  into  the  computer's  active  core  area  before 
execution  can  begin.  Loading  is  done  using  a command  file, 

BUS3.CMD,  which  contains  the  names  of  all  files  that  must  be  loaded. 
Section  5.6  contains  a listing  of  BUS3.CMD,  and  the  list  of  files 
to  be  loaded. 

Once  all  files  are  compiled,  the  following  command  is  typed 
on  the  user's  terminal: 

.EXECUTE  @ BUS3.CMD 

The  subroutines  are  then  loaded  along  with  the  main  program, 
and  execution  begins.  The  program  first  asks  the  user  to  specify 
one  of  three  typical  driving  cycles  which  determine  cruising  speed, 
deceleration  rate,  number  of  stops  per  mile,  route  length,  etc. 

A "1"  is  typed,  followed  by  a carriage  return  if  Cycle  A is  desired, 
a ”2"  for  Cycle  B,  or  a "3"  for  Cycle  C.  The  program  then  asks  for 
a code  number  to  indicate  which  type  of  data  is  desired.  A "1"  will 
cause  the  output  data  file  FOR03.DAT  to  be  generated.  (See  Section 
4.5.)  A ”2"  will  generate  the  plot  file  FOR22.DAT,  which  can 
later  be  used  by  the  plotting  program  to  produce  CALCOMP  plots  of 
system  parameters.  A "3"  will  generate  both  files,  a "4"  neither. 
The  code  number  is  typed,  followed  by  a carriage  return,  and  the 
program  proceeds  unassisted.  At  the  end  of  execution,  the  file 
FOR03.DAT  is  automatically  queued  for  printing  on  the  line  printer. 
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To  print  all  files  contained  in  the  command  file  BUSL3.CMD 
(Section  4.6),  the  following  command  is  given: 

.LIST  @ BUSL3.CMD  ■ 

3.3  PLOTTING  PROGRAM  CHART. FOR 

A plotting  program  has  been  developed  to  display  any  six 
time-varying  system  parameters  on  an  off-line  CALCOMP  pen  plotter. 
The  program  uses  the  standard  CALCOMP  subroutine  calls,  and  writes 
instructions  on  a magnetic  tape  which  are  later  interpreted  by  the 
plotting  device  to  generate  the  plots.  Figure  10  shows  a sample 
plot  for  Driving  Cycle  C.  The  peak  values  of  each  parameter  must 
be  determined  before  execution,  and  appropriate  y-axis  scales 
must  be  chosen.  These  scale  values  are  entered  during  execution  of 
the  program,  as  are  labels  for  each  of  the  six  strip-charts  and 
the  x-axis  length  in  inches.  Section  5.7  contains  a listing  of  the 
plotting  program,  which  has  the  file  name  CHART. FOR. 
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FIGURE  10.  CALCOMP  PLOT  OF  SYSTEM  PARAMETERS  - DRIVING  CYCLE 


H,  modeling  theory  and  equations 


This  section  details  the  modeling  theory  used  in  simulating 
the  propulsion  system's  components.  The  modeling  equations  are 
expressed  in  terms  of  FORTRAN  symbols  which  appear  later  in  the 
program  listing.  The  values  chosen  for  the  system  component  param- 
eters are  nominal  ones  selected  for  the  purpose  of  illustrating 
baseline  system  performance. 

4.1  MISSION  PROFILE  EQUATIONS 

The  vehicle's  mission  profile  is  defined  at  the  start  of  each 
run  by  selecting  one  of  three  typical  driving  cycles,  called 
cycles  A,  B,  and  C.  The  driving  cycles  are  completely  specified  by 
the  various  constants  given  in  Table  1.  The  actual  vehicle 
acceleration  (ACC)  is  determined  by  the  net  accelerating  thrust 
(TN)  supplied  by  the  propulsion  system.  During  acceleration,  a 
relay  in  the  cam  controller  maintains  the  motor  armature  current 
(AMP)  at  or  above  the  specified  setpoint  (SET) . This  is  done  by 
advancing  the  cam  controller  step  (ICC)  to  the  next  position  as 
the  current  falls  below  the  setpoint.  On  the  highest  step  (ICC=14) 
the  vehicle  is  operating  at  its  power  limit,  and  the  decreasing 
armature  current  is  allowed  to  fall  below  the  setpoint.  The  vehicle 
then  accelerates  at  a reduced  rate  up  to  the  cruise  speed  (VC) . 

The  jerk  (AJERK)  is  constrained  to  be  within  a specified  limit 
(AJERKR)  during  acceleration  and  deceleration.  During  cruise,  the 
jerk  limit  is  reduced  in  order  to  minimize  force  fluctuations  as 
the  cam  controller  maintains  cruise  speed  by  alternating  between 
motoring  and  braking  modes  of  operation. 

The  number  of  stops  per  mile  (NS)  specifies  the  distance 
between  stops  (SS) , since  equal  stop  spacing  is  assumed.  The 
dwell  time  at  each  stop  (DWELL)  is  the  time  allowed  for  passengers 
to  enter  or  exit  the  vehicle.  Additional  dwell  time  at  the  end  of 
the  route  (TAUMU)  can  be  added  to  simulate  a waiting  period  before 
the  next  transit  of  the  route  is  begun.  The  number  of  passengers 
(NP)  is  used  in  computing  the  vehicle's  gross  weight  (WT)  and 
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accelerating  mass  (AMASS) . The  route  length  (SR)  times  the  number 
of  transits  of  the  route  (NT)  gives  the  total  distance  traveled  by 
the  vehicle. 

The  deceleration  is  always  done  at  a constant  but  jerk-limited 
rate  (DECL) . For  low  braking  forces  at  higher  speeds,  regeneration 
to  the  trolley  wire  occurs  as  the  motor  is  switched  to  a generating 
configuration.  For  higher  braking  forces,  however,  the  regenerated 
power  is  diverted  from  the  trolley  wire  to  dynamic  resistors  on 
board  the  vehicle.  Any  additional  braking  power  to  meet  the 
required  deceleration  is  absorbed  by  friction  brakes.  Below  4.0 
mph,  friction  brakes  are  used  exclusively.  The  braking  is  begun  at 
a predetermined  distance  from  the  next  stop  so  that  the  vehicle 
comes  to  rest  at  the  proper  location.  Figure  3 shows  the  speed- 
time plot  for  one  complete  cycle  of  an  arbitrary  mission.  The  time 
segments  for  the  deceleration  are  determined  as  follows: 

1)  T1  = DECL/AJERKR  - Time  from  start  of  braking  to  start 

of  constant  deceleration  [sec] . 

2)  T2  = V/DECL  - Time  from  start  of  braking  to  completion  of 

constant  deceleration  [sec] . 

3)  T3  = T1  + T2  = Total  time  of  deceleration  [sec] . 

4)  T4  = T3  + DWELL  = Time  from  start  of  braking  to  end  of 

cycle  [sec] . 

5)  T4  = T5  + TAUMU  = Additional  time  added  between  routes. 
Where : 

DECL  = Constant  deceleration  rate  (mph/sec). 

? 

AJERKR  = Jerk  limit  (mph/sec  ) . 

V = Vehicle  speed  at  start  of  braking  (mph) . 

DWELL  = Dwell  time  at  each  stop  (sec) . 

TAUMU  = Dwell  time  between  routes  (sec)  . 

The  stopping  distance  (XB)  for  the  given  deceleration  rate 
(DECL)  and  vehicle  speed  (V)  is  computed  as  the  product  of  the 
total  time  of  deceleration  (T3)  times  the  average  vehicle  speed 
during  deceleration  (V/2.0).  That  is: 
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6)  XB  = (V/2.0)  * T3/3600.0  = stopping  distance  (mi). 


or : 

7)  XB  = (V/2.0)  * (DECL/AJERKR  + V/DECL)  / 3600.0 

At  each  program  iteration,  the  vehicle's  acceleration  is  com- 
puted as  the  net  accelerating  thrust  (TN)  divided  by  the  equivalent 
vehicle  mass  (AMASS).  In  computing  AMASS,  a factor  of  2300.0  lbs 
is  added  to  the  vehicle's  gross  weight  (WT)  to  account  for  rotary 
inertia  of  the  propulsion  system  components. 

8)  ACC  = TN/AMASS  = Vehicle  acceleration  (mph/sec) . 

9) AMASS  = (WT  + 2300.0)  / GEE  = Vehicle  accelerating  mass 

(g- lbs)  . 

10)  WT  = WTCRB  + (NP+1)  * WTP  = Vehicle  gross  weight  (lbf) . 
Where : 

TN  = Net  accelerating  thrust  (lbf) . 

GEE  = 21.927  = Acceleration  due  to  gravity  (mph/sec). 

WTCRB  = 22000.0  = Vehicle  curb  weight  (lbf). 

NP  = Number  of  passengers. 

WTP  = 150.0  = Weight  of  a typical  passenger  (lbf). 

The  average  vehicle  acceleration  (ABAR)  over  the  computing 
time  interval  (TAU)  is  then  determined  and  numerically  integrated 
to  determine  the  vehicle's  speed  (V)  during  that  time  interval. 

The  average  speed  (VBAR)  is  similarly  integrated  to  determine  the 
distance  traveled  by  the  vehicle  during  each  interval.  The  total 
time  since  the  start  of  the  run  (THETA)  is  the  sum  of  each  com- 
puting time  interval  (TAU).  During  acceleration,  cruise,  and 
deceleration,  TAU  equals  0.05  seconds;  but  during  dwell,  when  the 
system  parameters  vary  slowly,  TAU  is  increased  to  0.2  seconds  in 
order  to  reduce  the  total  computing  time.  The  total  distance 
traveled  by  the  vehicle  since  the  start  of  the  run  (S)  is  the  sum 
of  the  distances  traveled  during  each  iteration.  The  distance 
between  the  previous  stop  and  the  point  at  which  braking  (STAB) 
begins  equals  the  distance  between  stops  (SS)  minus  the  stopping 
distance  (XB) . When  the  vehicle  has  gone  the  distance  STAB, 
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braking  begins  at  the  predefined  deceleration  rate  (DECL) , with 
allowances  made  for  the  specified  jerk  limit  (AJERKR) . 

At  each  iteration,  the  drag  forces  acting  on  the  vehicle  are 
computed  according  to  the  following  equations: 

11)  TG  = WT  * SIN (AT AN  (GRADE) ) = Grade  drag  (lbf) . 

12)  TR  = CR  * WT  = Rolling  drag  (lbf). 

13)  TA  = 0.00258  * CD*AF* (V+HW) **2.0  = Aerodynamic  drag  (lbf) 

14)  TD  = TG+TR+TA  = Total  drag  force  (lbf). 

Where : 

WT  = Vehicle  gross  weight  (lbf) . 

GRADE  = Roadway  grade  (rad) . 

CR  = 0.0075  = Rolling  drag  coefficient. 

CD  = 0.66  = Aerodynamic  drag  coefficient. 

AF  = 69.5  = Vehicle  frontal  area  (ft^) . 

V = Vehicle  speed  (mph) . 

HW  = Encountered  headwind  (mph) . 

The  grade  data  is  stored  in  an  array,  G(25,2),  which  can 
contain  up  to  25  paired  sets  of  grade  values  and  the  corresponding 
distances  at  which  the  grades  are  encountered.  As  the  vehicle 
passes  each  distance  at  which  the  grade  changes,  a new  value  of 
GRADE  is  determined  from  the  next  indexed  value  in  the  array 
G(25,2).  In  this  simulation,  all  grade  values  are  zero. 

4.2  MOTOR  MODELING  EQUATIONS 

The  dc  traction  motor  has  a shunt-wound  field  that  is 
rheostatically  controlled.  The  motor  terminal  voltage  can  be 
varied  by  adding  or  removing  resistors  in  series  with  the  armature 
circuit.  A schematic  diagram  of  the  vehicle's  propulsion  system 
is  shown  in  Figure  11.  A 14-position  cam  controller  regulates  both 
series  and  shunt  resistance  values  in  the  motoring  and  regenerative 
braking  modes  of  operation.  During  acceleration,  a relay  maintains 
the  armature  current  (AMP)  at  or  above  the  specified  setpoint  (SET) 
This  is  done  by  advancing  the  cam  controller  to  the  next  position 
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MOTORING  AND  REGENERATIVE  BRAKING  MODES 


TROLLEY 


DYNAMIC  BRAKING  MODE 


FIGURE  11.  PROPULSION  SYSTEM  SCHEMATIC  DIAGRAMS 
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as  the  current  falls  below  the  setpoint.  Each  higher  cam  position 
removes  more  series  resistance  from  the  armature  circuit,  thus 
increasing  the  motor  terminal  voltage.  When  the  cam  position 
indicator  (ICC)  equals  10,  the  series  resistance  is  reduced  to 
zero.  The  field  weakening  mode  is  entered  at  position  11,  as  the 
field  resistance  is  increased  on  this  and  all  higher  cam  positions. 
On  the  highest  cam  position  (ICC=14)  the  power  limit  is  encountered, 
and  the  vehicle  accelerates  at  a reduced  rate  up  to  cruise  speed 
(VC) . Table  6 summarizes  the  values  of  series  resistance  (RS)  and 
shunt  resistance  (RSH)  for  all  14  positions  of  the  cam  controller. 

The  motor  speed  (RPM)  is  related  to  the  vehicle  speed  (V) 
by  the  following  equation: 

15)  RPM  = RATIO*TF*V=Motor  speed  (rpm) . 

Where : 

RATIO  = 11.59  = Reduction  gear  ratio 
TF  = 8.003  = Tire  factor  (rpm/mph) . 

V = Vehicle  speed  (mph) . 

The  field  current  (AMPF)  is  computed  as : 

16)  AMPF  = VOLT/RF  = Field  current  (amps) . 

Where : 

VOLT  = 550.0  = Trolley  wire  voltage  (volts). 

RF  = RFO  + RSH  = Field  resistance  (ohms) . 

RFO  = 83.0  = Constant  resistance  term  (ohms). 

RSH  = Variable  shunt  resistance  term  (ohms) . 

The  motor's  magnetization  curve  is  modeled  by  Equation  17, 
known  as  the  Froelich  Equation.  Equation  17  is  also  plotted  in 
Figure  12  for  the  14  cam  positions  summarized  in  Table  6. 

17)  EA  = EAO*  (AMP  + SIGMA1*AMPF)  * RPM 

(AMPO  + AMP  + SIGMA* AMPF)  1000.0 
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TABLE  6.  FROELICH  EQUATION  PARAMETERS 


Motoring  and  Regenerative  Braking 

ICC 

RS 

RSH 

EAO 

AMPO 

1 

3.197 

99999.9 

501.5 

252.9 

2 

3.197 

41.0 

562.1 

266.4 

3 

2.090 

41.0 

624 . 9 

266.4 

4 

1.451 

41 . 0 

624.9 

273.3 

5 

1.082 

41 . 0 

624.9 

273.3 

6 

0.812 

41.0 

624.9 

2 73.3 

7 

0.606 

41.0 

624.9 

273.3 

8 

0.358 

41 . 0 

624.9 

273.3 

9 

0.156 

41 . 0 

624 . 9 

273.3 

10 

0.0 

41 . 0 

624 . 9 

273.3 

11 

0.0 

103.0 

545.2 

212.1 

12 

0.0 

220  . 0 

515.1 

201 . 9 

13 

0.0 

520.0 

561 . 2 

266.4 

14 

0.0 

99999.9 

501.5 

252.9 

Rheostatic  Braking 

IREO 

RD 

RSH 

EAO 

AMPO 

1 

1.451 

103.0 

616.4 

336.0 

2 

1 .451 

41  . 0 

616.4 

336.0 

3 

1.451 

9.0 

616.4 

336 . 0 
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EA  (GENERATED  EMF) 


AMP/SIGMA1  + AMPF  ( EQ . FIELD  AHP) 


FIGURE  12.  DC  MOTOR  MAGNETIZATION  CURVES 
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Where : 


EA 
EAO 
AMP 
AMPF 
AMPO 
SIGMA 
SIGMA 
CAY PH  I 
RPM 


= Generated  emf  (volts). 

= Base  value  of  EA  (volts). 

= Armature  current  (amps)  . 

= Field  current  (amps). 

= Constant  current  term  (amps) . 

= 82.916  = Field/armature  turns  ratio. 

= SIGMA1  + CAYPHI  = Corrected  turns  ratio  term. 
= 15.0  = Flux  correction  term. 

= Motor  speed  (rpm)  . 


Also,  from  Figure  11: 

18)  EA  = VOLT  - VB -AMP*R  = Generated  emf  (volts). 
Where : 


VOLT  = 550.0  = Trolley  wire  voltage  (volts). 

VB  = 2.0  = Brush  voltage  drop  (volts). 

AMP  = Armature  current  (amps) . 

R = RA+RS  = Armature  circuit  resistance  (ohms). 

RA  = Motor  armature  resistance  (ohms) . 

RS  = Added  series  resistance  (ohms). 

Equating  Equations  17  and  18,  thereby  eliminating  the  variable 
EA,  and  solving  the  quadratic  expression  for  the  motor  armature 
current  (AMP)  yields  : 


19)  AMP 


-B+  /B**2.0  - 4.0  * C 

~ = Motor  armature  current 

Z-U  (amp). 


Where  : 


B = AMPO  + SIGMA* AMPF  - (VOLT-VB)/R 
+ EAO  * RPM/ (1000 . 0*R) 

C = -(VOLT-VB)  * (AMPO  + SIGMA*AMPF) /R 
+ EAO*RPM* SIGMA1 *AMPF/ (1000 . 0*R) 


Once  the  armature  current  (AMP)  is  determined,  the  terminal 
voltage  (VT) , generated  emf  (EA) , and  air-gap  power  (PI)  are  com- 
puted by  the  following  equations  : 
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20)  VT  = VOLT  - AMP*RS  = Motor  terminal  voltage  (volts) . 

21)  EA  = VT-VB-AMP*RA  = Generated  emf  (volts) . 

22)  PI  = EA*AMP/1000 . 0 = Air-gap  power  (kw) . 


NOTE:  For  regenerative  braking,  the  sign  of  the 

brush  voltage  drop  term  (VB)  must  be  reversed. 


term  (kw) . 

>M  = Bearing  loss  term  (kw) . 

P5  = 3.143  * (EA/550 . 0) **2 . 0 = Core  loss  term  (kw) 

i*RPM**3.0  = Friction  and  windage  los 
(kw)  . 


The 

motor 

loss  t 

23) 

P2  = 

3.  39  E 

24) 

P3  = 

7.156  : 

25) 

P4  = 

0 .001 

26) 

P5  = 

3.143 

27) 

P6  = 

1.77  E 

28) 

P 7 = 

P1-P2-: 

29) 

P8  = 

0.0099 

The  motor 

shaft 

power 

30) 

PS  = 

P7-P8  ’ 

The 

drive 

train 

equation : 

31) 

P9  = 

1.5  E- 

+ 

0.0652 

Where : 

RPM  = Motor  speed  (rpm) . 

RATIO  = 11.59  - Reduction  gear  ratio. 

PS  = Motor  shaft  power  (kw) . 

The  available  propulsion  power  (PW) , tractive  effort  (TE) , 
and  net  accelerating  force  (TN)  are  then  computed  as: 

32)  PW  = PS-P9  = Propulsion  power  (kw) . 

33)  TE  = PW/ (0 . 001989*V)  = Tractive  effort  (lbf) . 

34)  TN  = TE-TD  = Net  accelerating  force  (lbf) . 
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Where : 


PS  = Motor  shaft  power  (kw) . 

P9  = Drive  train  power  losses  (kw) . 

V = Vehicle  speed  (mph) . 

TD  = Total  drag  force  on  vehicle  (lbf) . 


For  rheostatic  braking,  the  magnetization  curve  is  modeled 
by  the  following  equation: 

35)  EA  = EAO*  ( ARH -SI GMA1 *AMPF ) * RPM 

( -AMPO+ARH- SIGMA* AMPF)  1000.0 


Where : 


EA 

EAO 

ARH 

AMPF 

AMPO 

SIGMA1 

SIGMA 

CAYPHI 

RPM 


Generated  emf  (volts) . 

Base  value  of  EA  (volts). 

Rheostatic  braking  current  (amps) . 

Field  current  (amps) . 

Constant  current  term  (amps) . 

82.916  = Field/armature  turns  ratio. 
SIGMA1+CAYPHI  = Corrected  turns  ratio  term. 
15.0  = Flux  correction  term. 

Motor  speed  (rpm) . 


Also,  from  Figure  11: 


36)  EA  = ARH* (RA+RD) +VB  = Generated  emf  in  rheostatic  braking 

mode  (volts) . 

Eliminating  EA  from  Equations  35  and  36  and  solving  the  re- 
sulting quadratic  for  the  rheostatic  braking  current  (ARH)  yields: 

37)  ARH  = -BB-  /BB**2.0  - 4.0*CC 

2 . 0 


Where : 

BB  = - (AMPO+S IGMA*AMPF -VB/ (RA+RD) 

+EAO*RPM/ (LOO  0 . 0* (RA+RD) ) ) 

CC  = - VB*  (AMPO+ SIGMA* AMPF)/ (RA+RD) 

+ EA0*RPM*SIGMA1*AMPF/ (1000 . 0* (RA+RD) ) 

Once  the  rheostatic  braking  current  (ARH)  is  determined,  the 
terminal  voltage  (VT) , generated  emf  (EA) , and  airgap  power  (PI) 
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are  computed  by  the  following  equations: 

58)  VT  = ARH*RD  = Terminal  voltage  in  rheostatic  braking  mode 

(volts) 

39)  EA  = VT+ARH*RA+VB  = Generated  emf  in  rheostatic  braking 

mode  (volts) . 

40)  PI  = -EA*ARH/1000 . 0 = Air-gap  power  in  rheostatic  braking 

mode  (kw) . 

Where : 

ARH  = Rheostatic  braking  current  (amps) . 

RD=1.451  = Rheostatic  braking  resistance  (ohms). 

RA=0.1267  = Armature  resistance  (ohms). 

VB=2.0  = Brush  voltage  drop  (volts). 

During  rheostatic  braking,  three  decreasing  values  of  field 
resistance  (RSH)  are  tried  in  order  to  obtain  the  required  braking 
thrust.  These  values  are  summarized  in  Table  6 along  with  the 
other  parameters  used  in  the  Froelich  Equation  (Equation  35).  If 
the  required  thrust  exceeds  that  available  with  the  rheostatic 
braking  system,  friction  brakes  supply  the  remainder  of  the  braking 
thrust  requirement.  The  motor  and  drive  train  loss  terms  apply  to 
both  motoring  and  braking  modes  of  operation.- 

4.5  PCU  MODELING  EQUATIONS 

The  power  conditioning  unit  (PCU)  loss  is  the  heating  power 
dissipated  in  the  series  resistance  (RS)  that  is  added  to  the 
armature  circuit  for  voltage  control.  Therefore: 

41)  PCUL  = RS*AMP**2 . 0/1000 . 0 = Power  conditioning  unit  loss 

(kw)  . 

The  motor  field  power  (PTMF)  is  computed  as  follows: 

42)  PTMF  = VOLT*AMPF/1000 . 0 = Traction  motor  field  power  (kw) . 

The  total  auxiliary  power  (PAUX)  includes  the  field  power 
(PTMF),  and  is  determined  by  the  following  equation: 

43)  PAUX  = PTMF+PECBLO+PAIR+ 

PLTG+PBC+PAC+PCBLO= 

Total  auxiliary  power  (kw) . 
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Where : 


RS 

AMP 

VOLT 

AMPF 

PECBLO 

PAIR 

PLTG 

PBC 

PAC 

PCBLO 


Armature  series  resistance  (ohm) . 

Motor  armature  current  (amp) . 

550.0  - Trolley  wire  voltage  (volt). 

Field  current  (amp) . 

0.6  = Environment  control  blower  power  (kw) . 

1.0  = Average  air  compressor  load  (kw) . 

3.6  = Lighting  power  (kw) . 

2.0  = Average  battery  charger  load  (kw) . 

6.0  = Air  conditioning  compressor  power  (kw) . 
0.6  = Air  conditioning  condensor  blower  (kw) . 


The  total  input  power  to  the  motor  (PTMA)  is: 


44)  PTMA  = VT*AMP/1000 . 0 = Motor  input  power  (kw) . 

The  total  input  power  to  the  vehicle  (PVEH)  is  then  computed 

as  : 


45)  PVEH  = PTMA  + PCUL  + PAUL  = Vehicle  input  power  (kw) . 
Where : 


VT 

AMP 

PCUL 

PAUX 


Motor  terminal  voltage  (volt) . 
Motor  armature  current  (amp) . 

Power  conditioning  unit  loss  (kw) . 
Total  auxiliary  power  (kw) . 
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5,  FORTRAN  SOURCE  LISTINGS  AND  DATA  FILES 
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5.1)  MAIN3 . FOR 


c***#  TROLLEY  BUS  PERFORMANCE  PROGRAM! 

COMMON  T A U 1 , lACr IUC, VC#N8,85»DWELL»TMJMUrNP#WT#SR#NT 

data  taui/i  ,o/»timx/o.2/#iac/o/ 

C****  SPECIFY  DRIVING  CYCLEf 
WRITE ( 5 , 2 1 4 ) 

READ  C 5 f 2 1 5 ) IDC 
C#***  SPECIFY  OUTPUT  FILES! 

WR ITE  C 5 , 2 1 8 ) 

READ  (5*215)  I0f  , 

IF (Clbr.EQ.il .OR. (Ior.EQ.35)  II  Si 
IF ( ( lor.EQ.2) .OR. ( IQF.EQ.3) ) 12*1 
WkITE  C 5 , 2 i 6 ) 
c**»*  INITIALIZATION | 

1*1 
11*2 
JMr  l 

CALL  PROF (V, ACC, SET, GRADE , HW , S , SN , TN , TD , TS , T4 , TAU , JM , 18W ) 

C****  INITIALIZE  CYCLE  DATAj 
35  ISW*1 
ACCsO.Ul 
VeC.Cl 

TAl'*TAUi/20.Q 
T I M e C . 0 
IWCsl 
SNsS^+SS 
ICCeO 

IF(Il.NE.l)  GO  TO  4u 
C***»  ntjrp'JT  P A GF  HEADING! 

WPITFC3,2l7) 

C***»  OUTPUT  DRIVING  CYCLE* 

IF(IDC.EQ.I)  WRITEC3.211) 

IF  f inc.EQ.23  WRITE(3*212) 

JFCI&C.EIO)  WRITE ( 3 , 2 1 3 ) 

C**»*  OUTPUT  AIR  CONDITIONING  STATUS! 

IF(IAC.EQ.t)  WRITE ( 3 * 20 1 ) 

IF(HC.EQ.O)  WRITE  f 3 # 202  ) 

C * ***  OUTPUT  CYCLE  DATA  j 

WRITE (3, 203)  I,II,SS,NP,WT, VC,DwELL,Hk 
C***»  OUTPUT  COLUMN  HEADINGS! 

WRITE ( 3 # 204 ) 

C**»*  DETERMINE  PROFILE  CONDITIONS! 

40  CALL  PR0F(V,ACC,SET,GRADF,Hw,s,SN,TN»TDrTs,T4,TAUf JM.ISO 
IF  fTS.GT.T4)  GO  TO  60  JEN!)  oF  CYCLE! 

C ****  COMPUTE  TRACTIVE  FORCE  AND  PROPULSION  POWER  REQUIREMENT! 

tp*tn*td 

PWso.OOl989*TP*V 

C****  COMPUTE  MOTOR  VOLTAGES,  CURRENTS,  AND  AVAILABLE  POWERj 

400  CALL  MnTnp(AMp(SFT,ARH,AuPF,VT,EA,V,VC,PS,PW,XF,TP,RS,IND,lCC) 


CONTINUED  ON  NEXT  PAGE 
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5.1)  MAIN . 3F0R  (CONT) 


C**»*  rr^pJT?-  MutOP  TERMINAL  POWER: 
eTMfc*VT*AMp/l  oOOw 

C ****  UPDATE  TRACTIVE  AND  ACCELERkTING  FORCES  j 
TFflND.EQ.l)  TP=T£ 

IFfl^D.FO. 1)  TN*TP-TD 

c***«  compute  friction  braking  forte: 
tfctd.Lt.tej  tm*tp-te 

IFf  I*'D.LE.T)  TMsC.O 

C * ***  COMPUTE  PCU  LOSS  and  AUXILIARY  pOWERj 
CALL  pCu(AMP,AMPF,VT,RS,PCULfPAuX,lAC) 
r***»  COMPUTE  total  input  POWER  TO  VEhICLEi 
PVFH*PTMA+PCUL+PAUX 
PRAk*(PVEH+PVEHi )/2.0 

pvehi«pveh 

c ****  compute  total  vehicle  energy : 

EVEH»EVEH+PBAR#TAU/3600. 5 
C***»  OUTPUT  TIME  HISTORY  OVER  CYCLE: 

IFf  (iWC.EQ.i).AND.(Ii  .Fa'.l))  wP.  ITE(3#206)  THE  r A , ACC  , V , S , GR  ADE 
* ICC , TND(TP,TEf TM,TSTD,PK, PS, AMP, A PH, VT,EA,PTMA , PVFH,EVEH 
IFf (TWC.FQ.n .AWT. f I2.EQ.1 ) ) WP I TE ( 2 ? * 2 1 9 ) THET A # V , TP , PW , PVEH 

« a*p,ea 

I wCsO 

C*.***  INCPEMENT  PRINT  CLOCK: 

tim«tim+tau 
IFf TIM.GE.TIMX3  IWCsl 
IFCTIM.GE.TIMX)  TIM.o.o 
c**»*  IDCREMEWT  PUN  TIUE  CLOCK: 

thet^*theta+tau 

IF (THtTA. GT. 100. 0)  STOP 
GO  TO  4o 

C * * **  END  OF  CYCLE: 
tr\  TS  = 0.0 
1 = 1 + 1 
11=1+1 

IFfS.GE.JM*SP-0.05)  GO  TO  65 
GO  TO  35 

?**+*  END  OF  ROUTE: 

6 K JMsjM+1 

IFfS.GE.NT*SR-0.n5)  GO  TO  7o 
GO  TO  35 

C***»  END  OF  MISSION: 

7p  ES»C*EVEH/S 

VAV0«3bOO.O*S/THETA 
T»IN«THETA/60. 0 
C**+*  WfTTF  PUN  SUMMARYj 

WPITE(3,20T)  S#TMIN,VAVE,EVEH,ESPC 
201  FORMATf///#  * AlP  CONDITIONING  On»,/) 

?02  FORMAT (///,'  AIR  CONDITIONING  OFF*,/) 
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5.1)  MAIN3 . FOR  (CONT) 


203 

( 

EORMATC//,*  From  STOP*, 14,* 
* DISTANCE  BETWEEN  STOPS! 

TO*, 14,//, 
SS«* 

» F7  , 3 , * 

MI',//, 

* 

* NUMRER  OF  PASSENGERS! 

Nps» 

,17,//, 

* 

* vehicle  Gross  weight! 

WTa* 

# F? , o , : 

IBS*,//, 

t 

; CRUISING  SPEED! 

vc«* 

# F7 , 1 , ' 

MPH*,//, 

* 

* DWELL  TIME! 

DWELL** 

, F7.ii' 

SEC*,//, 

t 

• ENCOUNTERED  HEADWIND! 

Hwb* 

» F7 . 1 , * 

MpH *,//) 

grade  icc  ind  t p 

PS  AMP  ARK  VT 


204  PORMATC///!  * THETA  ACC  V S 

# TE  TM  rN  TD  Pis1 

t EA  PTMA.  PVEH  FVEH',/) 

2 06  FOPMAT(  iX,F6. 1 ,2,Pfc.  1 .Ff.2,F6'.  3 , 2l  4 , 5r7 . 0 i eFtT.  1 , F6 . 2 ) 

207  FORMATC * RUN  S U*M ARY j * , / / , 

« 

I 
I 
A 


total  DISTANCE  TRAVELLED} 
* TIME  or  MISSION! 


* AVERAGE  vehicle  SPEEDi 
; TOTAL  energy  CONSUMED! 

* SPECIFIC  VEHICLE  ENERGY! 


S**,r7.2,*  Ml',//, 
TMIN«»,F7.2i  ' MJN  » , // , 
VAVE«SF7.2,  * MPH*,//, 
EVEHb  * , F7  » 2 , * Kw»HRS*,//, 
ESPC ■ * » F7 . 7 


KW-HRS/MI*,//) 


211  FCRMATC/,*  DRIVING  CYCLE  A * ) 

71?  FORMATC/,*  DRIVING  CYCLE  6*) 

213  rORM  AT ( / , * DRIVING  CYCLE  C') 

214  FORMAT  ( / , * ENTER  DRIVING  CYCLE!*,/,, 

I i rypg;  hj"  pop  a , "2"  FOR  R,  "3"  .FOR  Ct'f/) 

215  FDPMATCII) 

216  FORMAT ( / » * WORKING...*,//) 

? 1 7 FORMATC jHi , ///, * TROLLEY  BUS  PERFORMANCE  PROGRAM*,//) 
218  FORMAT  (/ , * ENTER  TYPE  OF  OUTPUT  DESIREDf',/, 

• i TY^L  " 1 " FHR  DATA  FILE  FORfi 3 . D A l , * , / , 

« * TYPE  *5"  FOR  PLOT  FILE  FOR22 . DAT , * , / , 

« * TYPE  "3"  FOP  BOTH  FILES,*,/, 

t ; TYPE  "4"  FOR  NEITHER  FILE!*,/) 

21°  EORMATC iX,7Flo. 3) 

STOP 


EnD 
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5.2)  PROF 3 .FOR 


SUBROUTINE  PROF (V, ACC , SET , GRADE , H w , S , SN , TN , TD , TS , T 4 , T AU , JM , ISW) 
DIMENSION  G ( 25  * 2 ) 

COMMON  T*U1  #TAC|IDC,  VC»N5#5SiDWELL#TAUMU#NP,WT.,SP#NT 
DATA  GEE/21 .9?7/,iJFRKR/3.5/,PHlR/0.0/#pHIW/0.0/»wS/O.U/# 

* AF/69 .5/,  CD/0.56/, CR/O. OC75/# WTCPB/22000. 0/#WTp/ 150.0/ 

IFfINI.EQ.l)  GO  TO  1000 
I n j b ! 

C**»*  INITIALIZATION t 

IFfinc.EQ.l)  GO  TO  1 
IFflOC.EQ.2)  GO  TO  7 
IFCIDC.EO.J)  GO  TO  3 
C***«  DRIVING  CYCLE  Ai 

1 ACCMS3.-5 

DECL* 1 .5 

VCB2O.0 

NS*5 

SS«1 .O/NS 
DwELLb 12.0 
TAUMUbO.O 
NP*20 
SRb b .0 
NT»4 

SETb2  QQ . 0 
GO  TO  4 

C***#  DRIVING  CYCLE  Bj 

2 ACCM* 3 , 5 
D F C L s 3 • 5 
VC  = 25 ,0 
NS»F 

SS» 1 .O/nS 
DWELLS  16 .0 
TAUMUbO.O 
NPsiC 
SRb6 . 0 
NT»4 

St Ts?50.0 
GO  TO  A 

C **♦*  DRIVING  CYCLF  Cr 

3 4CCM«3.0 
DECL*  3 • 5 
VCe2< .0 

N £ B 1 0 

SS*1 . 0 /NS 
DWELLbB . 0 
TAUMUbO.O 
NF*  1 5 
S R b 0 . 2 
WTsi 

S r"  T e ? 2 5 . 0 

C***«  COMPUTE  VEHICLE  GROSS  WEIGHT  AND  ACCELERATING  Mass: 


CONTINUED  ON  NEXT  PAGE 
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5.2)  PROF 3 . FOR  (CONT) 


i.  W'f *WTCRB+  ( MP+  1 ? #WTe 
4 KASSs(KT+73no. 05/GEE 

c***#  grade  data: 

G(l#l)*O.G 
G( \ ,2)«0.0 
GC?f 1 )s100.0 
G(?,2)S0.0 
*G»1 

AJFRKsAjFpKR 

RFTUPN 

C**»*  COMPUTE  vehicle  denamics: 

1000  ACG=TN/AMASS 

ABAP*(  ACC  + ACCn/2.0 
ACC  1 *ACC 
VeY+ARAr*TAU 
VBAP.e(V  + Vl  1/2.0 
VI  SV 

SsS  + VB4j,#IAU/3600.0 
C»***  DFTEPKInF  HEADWIND  AND  GPADEt 

Hw*.W5*cn5( o.oi7A5*(PHlh-PHIW) ) 

55  JGsKG+1 

TFfS.GE.GCJG#m  GO  TO  60 
GH  TD  as 
DO  KGsKG+1 
GU  TO  55 

65  GRADE*G C K G » 25 /I  00.0 
C ****  COMPUTE  VEH'ICLt  DRAG? 

TGaWT*SIR  CATAN  CGRADE) ) 

TR*CC •^T 

TAsu.0  0256*CD#AF#(V*Hw)**2.0*(V+.HiO/4BSCV  + Hfc  + o.001) 
TbsTR+TA+TG 
IFtV.LE.0.01)  TDso.O 
C#***  STOPPING  LOGIC  & SpEF-D  CONTROL: 

TNIbAJErKR*AMASS«TAU 
TPfISW.EQ.fy)  GO  TO  6*0 

XB«V*(DECL/AJE*KR+V/DECL) /7200 • o 

STAB»SN-XB 

IF  fS.GE.STAB)  GO  TO  640 

IFfV.GT.VC)  GO  T'n  HO 

IF(V.GE.0.97*VC)  GO  TD  100 

IPfV.GE. fVC-ACC**?.r / AJERKR/2  . <) ) ) GO  TO  95 

TUsTM+TKl 

IF  ( rfJ  .Gt.  AC  CM*  A HASS  ) TN=ACCM*AMaSS 
GO,  TO  120 

(•****  JP^FD  withIN  DFA^BANDj 
9b  tNbTN-TNI 

IF ctn.Lt.c.o)  TN*o.O 
GO  TO  120 

100  TFCTm.GT.o.O’  TNeT«. TNT/3.0 
irfTN.LT.o.05  TNsTN  + TM/3,c 


CONTINUED  ON  NEXT  PAGE 


43 


5.2)  PR0F3.F0R  (CONT) 


r,<!  to  120 

(!**♦#  OVEPSPtED; 

1 lb  TN«TN-TNI/3.0 

if  CTn.LT.-DECL*AMASS)  TN.-DECL*AM'ASS 
GO  TO  120 

c**»*  INITIATE  STOP* 

Mb  IF  CISW.EQ.5  ) GO  TO  650 
GO  TO  bbO 
65b  ACC  a 0 . 0 

TlaDECL/AJERKp 
TjaV  /DFCIj 
13*11+12 
T4*T  3+DwELL 

IF (S.GF.JM*SP -0.051  T4«T4+TAUMU 

ISW*0 

TS*0 . 0 

660  TS*TS+TAU 

IF (TS.GT.T1 ) GO  TO  6^0 
c * ***  onseti 

AJERKb-aJEPKR 
GO  TO  69? 

6 7 0 IFCTS.GT.T2D  GO  TO  680 

r#***  nr.Cfc.LERATlON » 

AJEpKsO.O 
ACCs-DEcL 
Gn  TO  691 

680  IFCTS.Gt.T3D  GO  to  fe9 0 
C***«  FLAIR: 

AJFrKsAJF RKr 
GO  TO  69? 

69c  TrCTS.GT.T4)  GO  TO  i?q 
(■***♦  OWELL; 

AjrpKso . 0 
ACC  = 0 . 0 
VsO. 001 

fN*u .0 

TDs 0 . 0 
SaSN 

TAPsTAU 1/5.0 
Gn  Tn  120 

69?  ACC=ACO+AJERK*TAU 
694  TN*AGC*AMASS 
l 2b  RETURN 
END 
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5.3)  M0T0R3 . FOR 


SUBROUTINE  MOTOR ( AMp, SET, ARH , AMpF , VT , E A , V , VC , PS , PW , TE , TP , RS , I ND , 

# ICC) 

DATA  PA/0. 1 26  V,RD/l.  451 /.RATI  0/11 .59  / , RFO/B  3 . 0 / , TF /8  . Or  3 / , 

« VOLT/550.0/,SIGMA/07.916/,CAYPHI/15.0/,SIGMA1/82.916/,VB/2.0/, 

I CAY2/3.3  9E-o4/,CAn/’>.i56E-oR/,CAY4/O.Oai/,CAY5/3.143/, 

« CAY6/1 .77E-10/,CAYR/0.00«9/,CAY9/1 . SE.ftb / , C AY  1 0 /6 . 5 25E-0 2 / , 

* PS 0/3. 197/ 

TFfV.GT.VC)  Go  TO  100 
IF CV.GT.C.01 ) GU  TO  Bo 

C****  VEHICLE  AT  REST  j 
INDbO 
TCCeO 
GO  TO  510 

50  IF CTP.GE.o.O)  GO  TO  300 
C*#*#  PRAKINGj 

IFfV.GE.4.0)  GO  TO  IOC 
C#***  FULL  FRICTION  BRAKING: 

T ND»4 

TCCeO 

C***#  TURN  OFF  MOTOR  IF  VEHICLE  IS  AT  PEST  Op  In  FRICTION  BRAKING  MODE 
510  AHPbO.O 
AFHsO . 0 
AMPF»0.0 
V T b o , 0 
EA*0,0 

ps*o  • 0 

TE*0.0 

PETURN 

100  IF(ICC.LT. 10)  INDb  3 

IFflND.EQ.3)  GO  TO  200 
DC  550  I NDX  = 2 , 3 , 1 
TNDbINDx 

IFCINd.EQ.3)  GO  TO  2°0 
c***#  rlgenefative  bpaking: 

DO  55c  ICCX*12, 10,-1 

ICCsICCx 

GO  Tn  10 

C***»  RhEOSTATlC  bPAKINGl 
200  nf-  550  IREQXb!  , 3,  i 
TPEO*IPEQX 
ICCsO 
GO  TO  3 

c***»  motoring: 

300  I N D s 1 

IF  f ICC . EO , 0 ) I CC * 1 
If fV.GT.0.97*VC)  ICC ■ 1 ft 
310  IFfICC.GE.10)  GO  TO  1ft 
IFflCC.GE.4)  GO  TO  4 
TF(ICC.FQ.I)  RSsPSO 
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5.3)  M0T0R3 . FOR  (CONT) 


4 


c*«»« 

to 


c***« 

3 


c***» 

20 


IFdCC.EO.t  ) 
TF{ ICC.EQ. t ) 
IFf ICC.EQ. t) 
IFCICC.EQ.2) 
IP(ICC.EQ,2> 
IFCICC.F0.2) 
tfcicc.eq,?) 
IF  (ICC.EQ, 1) 
IF (ICC.EQ, 3) 
IFriCC  .EG.  3) 
TFflCC.EQ.3) 


PS H =99999. 9 
EAo*591  ,5 
A^POx252.9 
PSsPS-i 
RSris  4 1 .0 
£A:=562. 1 
AMPr>*266  . 4 
RS»0 . 37*RSo 
RSH*  4 t .0 
EAo«5*2.  1 

AMPfi*2Pfe.  A 


GO  Tn  2 0 


IF  f ICC.EQ. 4J 
IFf ICC.EQ. 5) 
IF(ICC.EQ.6) 
TFf ICC.EQ.  9) 
IF(ICC.EQ.R) 
IF  c ICC »EQ, 9 ) 
RSMe4l.y 
Ea0s*24.9 
AI'Po*2'73.i 
GO  tQ  2 0 
MOTORING  AND 
IFriCC.fc-Q.  1 j) 


RS*o , 45  3 9*RS 0 

PS«n . 33R4*KSp 
kS*o , ?540*RSo 
RS*y . 1 P®b*RSO 
R5»0. t 1 20*RSu 
PS*O.f'4«8*RS0 


REGENERATIVE  BHAKJmGi 
RSHrli  . 


iracc.ro.  1 0)  EAo*f24.9 

IFflCC.FG. to)  AKpn*273.3 
IF(ICC.EU.II)  RSHsi o3,0 
IFflCC.EO.il1  EA0SS45.2 
IFjlCC.EQ.1.1  ) AMP  0*212. 1 
IFflCC.E0.12)  RSHB5/f).f) 
IFCICC.EQ.12)  FA^eSlS.l 
IFCICC.EQ.12)  AMP 0*201 .9 
TFCICC  .EQ.  t 3)  RSH  = 5 20 . 0 
IFCICC.EQ.13)  EA  0*5  b 1 ,2 
IF  C ICC. F. 0.13)  AMP0*256.4 
IF  (ICC.EQ.  M)  RSHs99<>99.9 
IFCICC.E0.14)  E A 0»5  0 1 .5 
TFriCC.EQ.14)  AMPo*252.9 


RS«0 . 0 

GO  TO  2 0 

PFi.EOSTAlIC  BRAKING  1 
IF(IRtQ.EO.i)  RSRb 1 0 3 . 0 
IF  clRtO.EQ.2)  P5H»4l.O 
IFCIRE0.EQ.3)  RSHbR.O 
EA0«M&.4 
AKP0*336.C 

COMPUTE  motor  speed  AND  FIELD  RESISTANCE  AND  CURRENT  j 

PPMbRATIO*tr#v 

RF*RFo+r5H 

AMPF*VOLT/RF 

IFCIND.E0.3)  GO  TO  30 
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5.3)  MOTORS . FOR  (TONT) 


C***«  MOTORING  and  regenerative  BRAKING} 

P bR  A +R  S 

E«AMPo+SlGMA*AMpr-CVOLT-VB) /PtEA0*RP^/( iOOO.O*P) 

CB-Cv°LT-VB)*(AMPo+SlGMA*AMpr) /r  + 

t EAO*RPy*SIGMAl*AMpF/( iOOC.O*R) 

C #•*■*•#  COMPUTE  armature  current  i 

AMpB ( -B  + SQPT  CB**2 . C-4 . 0«C  ) W2.0 
C»**«  COMPUTE  MOTOR  terminal  VOLTAGEi 
vt»volt-amp*rs 

C»**«  COMPUTE  AIR-GAP  EME  j 

EA*VT.AMP*RA-VB*AMP/AB?»C  amp + 0,001  ) 
c***«  COMPUTE  air-gap  POWER  j 
P1*EA*AmP/100o'.0 
ARHbO.O 

GO  to  4o 

Cfr**«  RHEOSTATIC  BRAKING! 

3n  Bb»-CAMpo+SIGMA*AMPF»VB/fRA#RDWEAO*RPM/(1006.0*fRATRD)  ) 5 

CC»-VB*(AMPo+SlGMA*AMpF)/CRATRD)t  . 

* EAO#PPM*SIGMAi*AMPF/(  J ft£>0 . 0*  C R A + RD  ) ) 

c **»*  compute  braking  current! 

ARHb(-B&-SQRTcPB**?.o-4.0*CC) ) / 2 , 0 

c**»*  COMPUTE  MOTOR  TERMINAL  VOLTAGE! 
vt*arh*rd 

C#***  COMPUTE  AIP-Gap  EMFj 
EA*VTtArH*PATVB 
C****  COMPUTE  AIR-GAP  POWER: 

Pl»-EA«ARh/l 000.0 
AM.peO.O 

c#***  COMPUTE  MOTOR  POWER  LOSSES! 

4ft  P2*C A Y2#RPM 

P3*CAY3#RPM#*2 .0 
PasC AY4*Rpm 

t>5aCAY5*(t  A/55  0.0  )**2.0 
PpeC AYB*RPm**3 . 0 
C7sP1 -P2-P3-P4-P5-P6 
®BsC AYR* ABS  C p7 3 

c *♦*«  COMPUTE  MOTOR  shaft  POWER; 

?S*P7-Py 

C**»*  COMPUTE  DRIVT  TRAIN  LOSSES! 

DQsCAY9*(Rt>M/pATI0)**2.0  + CAYlQ*ABSCPS) 

C***»  COMPUTE  opopULSION  pOWER  AND  TRACTIVE  EFFORT! 

BV sPS-PR 

Tfc«PW/O.OOt9B9/V 

IFCINL.EQ.1)  GO  TO  5B0  IMOTOrINGI 
IF  CIND.EQ.4)  RETURN  1 FULL  FRICTION  BpAK  I F’G  1 

IFfTE.LE.TP)  RETURN 
550  CONTINUE 

PWbO . 00 1 9R9*TP*V 
return 

5b0  IFf cV.Gt.o.97*VC) .AND. (TE.GE.TP)  3 RETURN 

TFCAMp.GE.SET)  RETURN 
IFflCC.LT. t 4)  ICCbICC+1 
IFflCC.LT.t4)  CO  TO  3 1 v 
return 
E N 0 
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5.4)  PCU3 . FOR 


SUBROUTINE  PCutAMp, AMpF, VT , RS , PCUL , P AUX , I AC) 

HAT  A PAC/6,ft/f  PCBLO/o.b/#PECBLO/0.fc/|PAlRM‘.0>»PtTG/3'.6/i 

* PHC/2.0/# VHLT/550.0/ 

C***«  COMPUTE  power  conditioning  UNIT  LOSSi 
PC  !'L*BS*  AMP  **2.0/1000.0 
C***«  COMplJlF.  MOTOR  FIELD  power  : 

PTMF*V0lT#AMPF/1 OOO.r 
C***»  COMPUTE  AUX IL I AH Y POWER J 

PAUX«PTMF+PECBLO+pAIRtPLTG+PBC 
iKf  IAC.F.Q.  1 ) PAUXsPATTX  + PAC*PCRI.'0 
RETURN 
END 
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TRHbl.tY  A II 5 PFRFOFMa'JCp:  PHPtjPA'' 


5.5)  OUTPUT  DATA  FILE  - FOR0.3DAT 


PARTIAL  LISTING 


u 


X 

UJ 

c 

o 

© 

o 

fv 

© 

© 

o 

r 

o 

* 

-y 

3 

t~ 

o 

a o 

O 3> 
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- 

fN 

fN 

V 

<■ 

P- 

a- 

o 

c 

c 

rv 

rv 

> 
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o 

© 

© 

O 
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° 

o 

o 

O 

•D 
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o 

O 

c. 

o 

O 

O 

X 

O 

o 

O 

o 

o 

o 

O 

3 

rv 

X 

m 

c 

fN 

— 0* 

c 

CN 

fN 

a 

a 

fN 
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fN 
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> 

r— 

fN 

r- 

c 

o 

c 

c 

fN 

© 

T~  T~ 

X 
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W‘ 
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a. 

CD 

x 
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m 
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© 
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a 

fN 
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o 

If 
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•f 

C 

c 

•- 

r~ 
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c 

c 

a o 

a 

if 

a 
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fN 

T7 

ir 

V 
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Cr 
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IN 
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5.5)  OUTPUT  DATA  FILE 

RUN  SUMMARY: 

TOTAL  DISTANCE  TRAVELLED: 
TIME  OF  MISSION: 

AVERAGE  VEHICLE  SPEED: 
TOTAL  ENERGY  CONSUMED: 
SPECIFIC  VEHICLE  ENERGY: 


PARTIAL  LISTING  (CONT) 

S = 0.30  MI 
TMIN  = 1.59  MIN 
VAVE  = 11.34  MPH 
EVEH  = 1.30  KW-HRS 
ESPC  =4.33  KW-HRS /MI 


FORO. 3DAT 


5.6)  COMMAND  FILES  - BUS3.CMD  AND  BUSL3.CMD 


M4IN3.REL 
PROFS. REL 
PCU3.REL 
MOTORS . REL 


MAIN3 . FOR 
,PROF3 .FOR 
, PCU3 . FOR 
, MOTORS. FOR 
, F0R03.DAT 
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APPENDIX 


GLOSSARY  OF  PROGRAM  CONSTANTS  AND  VARIABLES 
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SYMBOL 

DESCRIPTION 

UNITS 

ABAR 

Average  acceleration  over  TAU. 

mph/sec 

ACC 

Vehicle  acceleration  rate. 

mph/sec 

ACC1 

Previous  value  of  ACC. 

mph/ sec 

AC  CM 

Thrust  limited  acceleration. 

mph/ sec 

AF 

Vehicle  frontal  area. 

ft2 

AJERK 

Jerk  rate. 

mph/ sec“ 

AJERKR 

Maximum  jerk  rate. 

mph/sec^ 

AMASS 

Vehicle  accelerating  mass. 

g-  lb 

AMP 

Motor  armature  current. 

amp 

AMPF 

Motor  field  current. 

amp 

AMPO 

Current  term  in  Froelich  Equation. 

amp 

ARH 

Rheostatic  braking  current. 

amp 

B,  BB 

Quadratic  equation  factors. 

-- 

C,  CC 

Quadratic  equation  factors. 

-- 

CAYPHI 

Flux  correction  term  in  Froelich  Equation. 

-- 

CD 

Aerodynamic  drag  coefficient. 

CR 

Rolling  drag  coefficient. 

C2-C10 

Motor  loss  term  coefficients. 

DECL 

Vehicle  deceleration  rate. 

mph/sec 

DWELL 

Dwell  time  at  each  stop. 

sec 

EA 

Generated  emf. 

volt 

EAO 

Base  value  of  EA. 

volt 

ESPC 

Specific  vehicle  energy. 

kw-hr/rai 

EVEH 

Vehicle  input  energy. 

kw-hr 

G (25  ,2) 

Grade  data  array. 

GEE 

Acceleration  due  to  gravity. 

mph/ sec 
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SYMBOL 

DESCRIPTION 

UNITS 

GRADE 

Roadway  grade. 

rad 

HW 

Encountered  headwind. 

mph 

I, II 

Stop  indices. 

— 

I AC 

Air  conditioning  status  code. 

— 

ICC 

Cam  controller  step  index. 

— 

IDC 

Driving  cycle  code. 

— 

IND 

Operating  mode  index. 

— 

IND  = 0 

Vehicle  at  rest. 

— 

IND  = 1 

Motoring  mode. 

— 

IND  = 2 

Regenerative  braking  mode. 

— 

IND  = 3 

Rheostatic  braking  mode. 

— 

IND  = 4 

Full  friction  braking  mode. 

— 

INI 

Initialization  bypass  code. 

— 

IOF 

Output  file  selection  code. 

— 

IREQ 

Rheostatic  braking  control  code. 

— 

ISW 

Stopping  logic  control  code. 

— 

IWC 

Output  write  command  code. 

— 

11,12 

Output  control  codes. 

— 

JG 

Grade  index. 

— 

JM 

Route  length  mile  index. 

— 

KG 

Grade  index. 

— 

NP 

Number  of  passengers. 

— 

NS 

Number  of  stops  per  mile. 

— 

NT 

Number  of  transits  of  route. 

— 

PAC 

Air  conditioning  compressor  power. 

kw 

PAIR 

Average  air  compressor  load. 

kw 

55 


SYMBOL 

DESCRIPTION 

UNITS 

PAUX 

Auxiliary  power  load. 

kw 

PBAR 

Average  vehicle  power  over  TAU. 

kw 

PBC 

Average  battery  charger  load. 

kw 

PCBLO 

Air  conditioning  condensor  blower. 

kw 

PCUL 

Power  conditioning  unit  loss. 

kw 

PECBLO 

Environmental  control  blower. 

kw 

PHIR 

Route  direction. 

deg 

PHIW 

Wind  direction. 

deg 

PLTG 

Lighting  power. 

kw 

PS 

Traction  motor  shaft  power. 

kw 

PTMA 

Traction  motor  armature  power. 

kw 

PTMF 

Traction  motor  field  power. 

kw 

PVEH 

Vehicle  input  power. 

kw 

PVEH1 

Previous  value  of  PVEH. 

kw 

PW 

Propulsion  power. 

kw 

PI 

Air  gap  power. 

kw 

P2-P8 

Traction  motor  power  loss  terms. 

kw 

P9 

Drive  train  losses. 

kw 

R 

Armature  circuit  resistance. 

ohm 

RA 

Armature  resistance. 

ohm 

RATIO 

Reudction  gear  ratio. 

RD 

Dynamic  braking  resistance. 

ohm 

RF 

Field  resistance. 

ohm 

RFO 

Minimum  value  of  RF. 

ohm 

RPM 

Traction  motor  speed. 

rpm 

RS 

Armature  circuit  series  resistance. 

ohm 
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SYMBOL 

DESCRIPTION 

UNITS 

RSO 

Series  resistance  factor. 

ohm 

RSH 

Shunt  field  resistance. 

ohm 

S 

Distance  traveled  since  start  of  run. 

mi 

SET 

Armature  current  setpoint. 

amp 

SIGMA 

Corrected  turns  ratio. 

SIGMA1 

Field/Armature  turns  ratio. 

SN 

Location  of  next  stop. 

mi 

SR 

Route  length. 

mi 

SS 

Distance  between  stops. 

mi 

STAB 

Location  of  brake  initiation  point. 

mi 

TA 

Aerodynamic  drag. 

lbf 

TAU 

Computing  interval. 

sec 

TAU1 

Factor  for  computing  TAU. 

sec 

TAUMU 

Time  added  at  end  of  route. 

sec 

TD 

Total  drag  force. 

lbf 

TE 

Tractive  effort  of  propulsion  system. 

lbf 

TF 

Tire  factor. 

rpm/mph 

TG 

Grade  drag. 

lbf 

THETA 

Time  from  start  of  run. 

sec 

TIM 

Output  time  clock. 

sec 

TIMX 

Output  interval. 

sec 

TM 

Mechanical  (friction)  braking  force. 

lbf 

TN 

Net  accelerating  force. 

lbf 

TNI 

Jerk  limited  thrust  increment. 

lbf 

TP 

Tractive  force. 

lbf 

TR 

Rolling  drag. 

lbf 
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SYMBOL 

DESCRIPTION 

UNITS 

TS 

Stopping 

clock . 

sec 

T1 

Time  to 

completion  of 

onset . 

sec 

T2 

Time  to 

completion  of 

constant  deceleration. 

sec 

T3 

Time  to 

completion  of 

flair . 

sec 

T4 

Time  to 

completion  of 

dwell . 

sec 

V 

Vehicle 

speed . 

mph 

VAVE 

Average 

vehicle  speed 

over  mission. 

mph 

VB 

Brush  voltage  drop. 

volt 

VBAR 

Average 

vehicle  speed 

over  TAU. 

mph 

VC 

Cruise  speed. 

mph 

VOLT 

Trolley 

wire  voltage. 

volt 

VT 

Traction 

motor  terminal  voltage. 

volt 

WS 

Wind  speed. 

mph 

WT 

Vehicle 

gross  weight. 

Ibf 

WTCRB 

Vehicle 

curb  weight. 

lbf 

WTP 

Weight  of  typical  passenger. 

lbf 

XB 

Stopping 

distance . 

mi 

110  Copies 
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